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A simple mathematical homogeneous model of competition is used to de- 
scribe Neanderthal extinction in Europe. It considers two interacting species, 
Neanderthals and Early Modern Men, in the same ecological niche. Using 
paleontological data we claim that the parameter of similarity, between both 
species, fluctuates between 0.992 and 0.997. An extension of the model in- 
cluding migration (diffusion) is also discussed nevertheless, extinction of Ne- 
anderthal seems unavoidable. Numerical analysis of travelling wave solutions 
(fronts) confirms the extinction. The wave-front-velocity is estimated from 
linear analysis and numerical simulations confirm this estimation. We con- 
jecture a mathematical formulation for the principle of exclusion between 
competitive interacting species (Gause). 
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Mathematical models for interaction and coexistence between species, 
yield non-linear equations which contain a number of rate constants. These 
constants are usually determined empirically in controlled experiments, or 
by observations in natural environments [1-5]. 

In this work we consider a system with competitive exclusion. We study 
the competition between Neanderthals and men with early modern structure. 
So, explicitly, we assume that Neanderthals did not evolve into Early Modern 
Men. 



Neanderthals were very stable in Europe for more than 60.000 years nev- 
ertheless, 40.000 years before our epoch, they were replaced by Early Modern 
Men. This massive extinction was completed between 5.000 and 10.000 years 
depending on the region. It is not clear how modern men appear in Europe, 
but some evidence of coexistence in the Levant (see for instance [6] and [7]), 
suggests a diffusive process from that region. Moreover, whether the interac- 
tion was direct (war) or not (competition) remains unclear nevertheless, we 
assume the last possibility. So, we have two species competing for the same 
ecological niche composed of limited nutrient and territory. The principle of 
exclusion (Gause) [5,8] can be applied and then both species cannot coexist. 



The following simplified non-linear mathematical model could describe 
this biological interaction. Let N be the number of individual, at time t, 
with Neanderthal characters, and C the one related to Early Modern Men. 
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Consider the pair of time-evolution equations 



dN 
~dt 



N(F-(3), 



(1) 



dC 
~dt 



C(F-sf3). 



(2) 



In this way, we assume the same rate of birth F(N, C) for both species, and 
a tiny difference in the constant rate of mortality (3. Namely, we assume 
that the parameter of similarity s varies between < s < 1, where s = 1 
means complete similarity. Thus, species C is better adapted to survive. The 
limited nutrient reservoir and territory is modeled by the dependence of F 
on the variables iV and C (see for instance[l]) : 



where the growth rate a and the interaction parameter 8 are positive num- 
bers, and we assume a > (5 > 0. 

Equations (1-3) define a Malthusian-birth-death process with interaction, 
and obviously is a crude model. For instance, it does not consider the dif- 
fusion process due to the alleged migration of modern hominid. Moreover, 
one might expect the rate constant a, f3 and 5 to be affected by changes of 
ecological nature (temperature, humidity, random fluctuations, amount of 
resources, and others). In fact, usually in population dynamics one deals 
with systems subjected to random environment [2,3] where, sometimes, the 
impact may be drastic (noise-induced transition). However, equations (1- 
3) are in accord with the principle of exclusion. Predictions related to the 
degree of similarity s, between both species, will be conjectured. 



F = a - 5(N + C) 



(3) 
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The set (1-3) can be solved partially. The usual linear stability analysis 
tell us that the point (N,C) = (0,0) is a unstable node; (N, C) = (^,0) 
is a saddle point (i.e. unstable); and finally (N,C) = (0, a ~ s ) is a stable 
node. All this in accord with Gause's principle (figure 1). Moreover, a direct 
integration gives : 



where A is a constant of integration. As long as < s < 1 (already assumed) 
the species N disappear, and the time of extinction r can be related to the 
parameters by 



Using the above relationship, paleontological data for the extinction time r 
(i.e. 5000 < r < 10000, years), and the life-time for individual (30 < 1//3 < 
40, years), then s fluctuates between 0.992 < s < 0,997. It is instructive to 
compare this result, for instance, with the parameter of similarity between 
man and chimpanzee [9] where s' ~ 0.975. This last parameter is related to 
the sequence of nucleotides in DNA, and then not necessarily connected to 
s. On the other hand, following Reef. [10], Neanderthal industry (silex-knife) 
did require about 111 percussions (4 stages) against 251 (9 stages) for Cro- 
Magnon. So the comparison of the number of percussion (for stage) gives 
s" — ^V^jp ~ 0.995. A number curiously close of our parameter s. 

Turning to the two species model (1-3), a more realistic case requires 
migration. This can be carried-up by adding a diffusive term onto (2), and 
neglecting the mobility of N in a first approach. In a more quantitative form, 
for instance, consider a migration term like to this one considered originally 
by Volterra [1,2], i.e. adding up a positive constant m onto (2) : 




■P(l-s)t 



(4) 



dC 
~dt 



C(F - s/3) + m, 



(2') 
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the linear analysis of (1,2') and (3) shows that the only stable point is 
(N,C) ~ (0, + — ), corresponding to the extinction for the species N. 
Nevertheless, a modification on the time of extinction r (5) is expected in this 
case. So, a constant migration term does not stop Neanderthal extinction. 



In a more realistic approach, we can consider a diffusive term like to 
Dd xx C added to (2) and look for travelling wavefront solution. Namely 
we consider solutions like to N(x — vt) and C(x — vt) in the population 
variables where v is the velocity of propagation. In this case using the variable 
z = x — vt , the evolution equations become, 

-v^ = N(a-(3-5(N + C)), (6) 

-v^ = C(a-sP-5(N + C)) + D^. (7) 

The linear stability analysis tell us that the point (N,C) = (^^ ,0) is a 
saddle point with one stable manifold. Nevertheless, the condition 



v 2 > AD (1 - s) (3 (8) 

is necessary because any physical solution requires C > 0. On the other 
hand the point (N, C) = (0, ^-^) is a unstable node and, finally, the point 
(N,C) = (0,0) is always stable (figure 2). For this last point, the condition 



v 2 > AD (a — sf3) (9) 

must be imposed (N, C > 0). Moreover, remark the invariance under velocity 
inversion (y — > —v) and coordinate inversion (z — > —z) in (6,7). So, for any 
solution with velocity v we can found one other with velocity —v. 
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From the above discussion, there is the possibility of travelling wave so- 
lutions connecting the point (^^,0) and (0,^-^). Numerical calculations 
confirm this possibility. Figure 3 shows two front-expansion for species C and 
the backward motion for iV is superposed. The stability of this wavefront 
solution was tested numerically using different extended initials conditions 
for N. So, numerical solutions confirm the extinction. Remark, like to the 
Fisher equation [5], we can use the inequality (9), more stronger than (8), to 
obtain a first estimation of the wavefront velocity, i.e. v ~ 
Numerical calculations confirm this estimation. Figure 4 shows the variation 
of the velocity wave front for different values of \[T). It was carried up as- 
suming / C(x, t)dx oc vt after the transient producing the stable front (figure 
3). 



Finally, the deterministic model (1-3) with interaction is simple and sim- 
ilar to those proposed in [1]. Nevertheless, as long as we assume the validity 
of the exclusion principle (Gause), other deterministic models [2-4] cannot 
gives us very different results. Namely, the exclusion principle guaranties the 
extinction of species N. In fact, from equations (1-3) we have the inequality 

ldC I ON 

clft > nW (10) 

which tell us that no-equilibrium points exist like to (N , Co), where iVo 7^ 
and Co 7^ 0. Moreover, the above inequality (10) guaranties the instability 
of equilibrium points like to (iVo,0). In this way, the above inequality can 
be interpreted as a mathematical formulation of the exclusion principle, and 
then we can expect its validity for other deterministic models describing the 
interaction between N and C. 
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FIGURE CAPTIONS 



Figure 1: A numerical sketch of competition between Neanderthal (N) 
and Early Modern Men (C). A simple mathematical model is described by 
equations (1-3). Species C fills gradually the ecological niche of N (arbitrary 
unities). 



Figure 2: A linear stability analysis of equation (6-7) with a diffusive term. 
(N, C) = (2=2, 0) is a saddle point and (N, C) = (0, 2r«2) is a unstable node. 
In this way, there is a possible solution connecting these points. The point 
(0, 0) is stable. 

Figure 3: A typical numerical simulation of the solution of (1-3) with a 
diffusive term Dd^C add to (2). The central pick corresponds to the expand- 
ing front for C. The other two, are related to Neanderthal backward motion 
(0 < X < 1000). The stability of these wave- fronts have been tested numer- 
ically using different extended initial conditions for N (arbitrary unities). 

Figure 4: A numerical calculation of the wave front velocity v for different 
diffusion coefficients (y/~D). The assumption, to obtain this graphic, is the 
number of total individual C grows like to vt after the transient (i.e. t — > oo). 
This calculation confirms our estimation for the wave front velocity. 
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